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Abstract: Otindag Sandy Land in China is an important ecological barrier to Beijing; the changes in its 
ecological quality are major concerns for sustainable development and planning of this area. Based on 
principal component analysis and path analysis, we first generated a modified remote sensing ecological 
index (MRSEI) coupled with satellite and ground observational data during 2001—2020 that integrated four 
local indicators (greenness, wetness, and heatness that reflect vegetation status, water, and heat conditions, 
respectively, as well as soil erosion). Then, we assessed the ecological quality in Otindag Sandy Land during 
2001-2020 based on the MRSEI at different time scales (i.e., the whole year, growing season, and 
non-growing season). MRSEI generally increased with an upward rate of 0.006/a during 2001-2020, with 
clear seasonal and spatial variations. Ecological quality was significantly improved in most regions of 
Otindag Sandy Land but degraded in the southern part. Regions with ecological degradation expanded to 
18.64% of the total area in the non-growing season. The area with the worst grade of MRSEI shrunk by 
15.83% of the total area from 2001 to 2020, while the area with the best grade of MRSEI increased by 
9.77% of the total area. The temporal heterogeneity of ecological conditions indicated that the 
improvement process of ecological quality in the growing season may be interrupted or deteriorated in the 
following non-growing season. The implementation of ecological restoration measures in Otindag Sandy 
Land should not ignore the seasonal characteristics and spatial heterogeneity of local ecological quality. The 
results can explore the effectiveness of ecological restoration and provide scientific guides on sustainable 
development measures for drylands. 
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1 Introduction 


Otindag Sandy Land in China presents a crucial ecological barrier and is a typical and fragile 
dryland ecosystem adjacent to Beijing (the capital of China) that is suffering from structural 
damage to parafunction, leading to land degradation and biodiversity loss (Wang et al., 2017; Wu 
et al., 2019). To slow down the process of land degradation and improve the ecological 
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environment, China has implemented a series of ecological restoration projects, such as the 
Three-North Shelter Forest Program and the Beijing-Tianjin sand source control project (Wu et 
al., 2013; Sun et al., 2019; Wen et al., 2020). Nevertheless, there is a continued debate on the 
effectiveness of these ecological restoration projects that have been implemented for 20 years 
(Yuan et al., 2016; Ma et al., 2017; Sun et al., 2019; Cong et al., 2022). A quantitative assessment 
of spatiotemporal dynamics in ecological quality is essential for the sustainable development of 
Otindag Sandy Land. 

With the rapid development of remote sensing and GIS technology, several remote 
sensing-based indices had been proposed to present the characteristics of ecological quality (Xu, 
2013; Xu et al., 2018; Wei et al., 2020). For instance, Wei et al. (2020) calculated the 
environmental vulnerability distance index by considering hydro-meteorological, socio-economic, 
soil-biological, and topographical factors to assess the ecological quality of Shiyang River Basin in 
China. Xu et al. (2018) developed a remote sensing ecological index (RSEI) based on the principal 
component analysis (PCA) to avoid the artificial influences on weights and yield more objective 
and reasonable; thus, this index has been widely used in the regional ecological quality assessment. 
However, researchers also realized that RSEI might have problems in ecological quality 
assessment in areas where urbanization may not be the main factor influencing ecological quality 
(Wang et al., 2020; Zhu et al., 2021; Yao et al., 2022; Ye and Kuang, 2022; Zhao et al., 2022). 
Therefore, the traditional RSEI has been improved and expanded in recent years. For example, 
Wang et al. (2020) developed an arid remote sensing ecological index (ARSEI) by replacing the 
dryness indicator in RSEI with the indicators of salinity and land degradation. Yao et al. (2022) 
modified RSEI with bare soil index and salinity index to reflect the macro ecological quality in 
arid areas. For applications in rocky desertified areas, Ye and Kuang (2022) constructed a modified 
remote sensing ecological index (MRSEI) by combining the degree of rocky desertification. An 
improved composite remote sensing ecological index (CRSEI) was developed by incorporating 
population density for studying the ecological quality of Yellow River Delta, China (Zhao et al., 
2022). In addition, to adapt to the local contradictions, researchers have developed new RSEI that 
can accurately characterize the changes in regional ecological quality, such as the remote sensing 
ecological index of local adaptation improvement (RSEILA) (Zhu et al., 2021). 

Ecological quality reflects the comprehensive characteristics of the elements, structures and 
functions of an ecosystem (Wu et al., 2019). Dryland is a special ecosystem type with sparse 
precipitation and high evaporation (Wu et al., 2019). In dryland ecosystem, vegetation is an 
important guarantee of regional ecosystem stability and greening is a major measure to mitigate 
desertification (Mamat et al., 2018; Chen et al., 2022). Meanwhile, water and heat are the principal 
natural conditions for vegetation growth and ecosystem development (Wang et al., 2021). Dryland 
ecosystem is the most fragile ecosystem, which is extremely vulnerable to human activities and 
climate change (Chen et al., 2022; Cong et al., 2022). Thus, soil erosion is a major reason for soil 
degradation and desertification in drylands that may augment when soil erosion proceeds faster than 
pedogenesis (IPCC, 2019). Soil erosion also reduces the capacity of water and soil conservation, 
affects regional water reserves, and threatens ecological quality seriously (Masroor et al., 2022). In 
this study, with a focus on comprehensive ecological quality assessment based on the regional 
characteristics of drylands, we combined the vegetation status, water and heat conditions, and soil 
erosion into an integrated index to assess the ecological quality in Otindag Sandy Land. 

As aforementioned, remote sensing image processing technology is flourishing; ecological 
indices based on remote sensing technology have been increasingly used in ecological quality 
assessment. However, adverse atmospheric states, sensor failures, and revisit periods still 
inevitably restrict the continuity of remote sensing data and hinder the wide application of remote 
sensing-based ecological indices (Qureshi et al., 2020; Yang et al., 2022a). Multiple data 
reconstruction methods aiming at improving the image quality, such as the maximum value 
composite, local filtering, and function fitting (Priyadarshi et al., 2018; Yang et al., 2022a, b; 
Zheng et al., 2022), are sound candidates, albeit these methods are highly dependent on the 
reliability of remote sensing data. Therefore, based on the PCA and path analysis of four 
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representative indicators (greenness, wetness, heatness, and soil erosion) of ecological quality in 
Otindag Sandy Land, this study attempted to construct a modified RSEI (i.e., MRSEI) coupled 
with satellite and ground observational data. Specifically, a long-time series MRSEI data 
(2001-2020) at different time scales (i.e., the whole year, growing season, and non-growing 
season) were calculated and used to assess ecological quality and its spatiotemporal variations in 
Otindag Sandy Land. The results can explore the effectiveness of ecological restoration and 
provide scientific guides on sustainable development measures in Otindag Sandy Land as well as 
other similar dryland regions in the world. 


2 Materials and methods 


2.1 Study area 


Otindag Sandy Land (41°46’-45°69'N, 111°55’-118°38’E) is located in the central part of Inner 
Mongolia Autonomous Region, China, and is one of the major sandy areas in the world (John et 
al., 2009; Ma et al., 2017). The landform of this area is relatively undulating, with a roughly 
higher terrain in the southeast and a lower terrain in the northwest, covering an area of 
approximately 0.17x10° km? (Liu and Wang, 2010). Otindag Sandy Land belongs to the arid and 
semi-arid continental climate zone, with the annual average temperature of 0.9°C—5.5°C. Affected 
by the monsoon, precipitation is gradually becoming less from the southeast to the northwest, 
with annual precipitation decreasing from 350-400 mm to 100-200 mm (Lou et al., 2019). 
Grassland is the major land cover type distributed across the area, accounting for 81.37% of the 
total area in 2020. Forests are patchy (5.13%) in the east while croplands are found in the south 
and east, occupying 1.92% of the total area (Fig. la). From 2001 to 2020, the land cover 
transformation occurred in 17.97% of the total area, of which 11.35% was transformed from 
grassland to cropland, 27.03% from bare land to grassland, and 11.18% from cropland to 
grassland, while 6.18% represented afforestation (Fig. 1b). 
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Fig. 1 Spatial distribution of different land cover types in Otindag Sandy Land in 2020 (a) and land cover 
transformation from 2001 to 2020 (b). The 1-km gridded land use cover datasets in 2000 and 2020 were derived 
from the Resource and Environment Science and Data Center (http://www.resdc.cn). 


2.2 Data sources and processing 


Monthly Normalized Difference Vegetation Index (NDVI) products (MOD13A3) from the 
Moderate-resolution Imaging Spectroradiometer (MODIS) at a spatial resolution of 1 km (from 
January 2001 to December 2020) were collected from the National Aeronautics and Space 
Administration (NASA) (https://ladsweb.modaps.eosdis.nasa.gov). NDVI data were used to 
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calculate greenness index (Zg). 

The daily meteorological data (2001-2020) of 103 weather stations, including the average air 
temperature (°C), maximum air temperature (°C), minimum air temperature (°C), relative 
humidity (%), sunshine duration (h), wind speed (m/s), and precipitation (mm) were from the 
National Meteorological Science Data Center (http://data.cma.cn). These data have undergone 
quality control (e.g., outlier rejection) to ensure their authenticity and reasonableness. We 
interpolated the daily data after quality control based on the Australian National University Spline 
(ANUSPLIN) software (Hutchinson and Xu, 2013), and extracted the final 1-km spatial 
meteorological datasets of Otindag Sandy Land for calculating wetness index (/,), heatness index 
(In), and soil erosion index (/;). ANUSPLIN has the advantages of solid theory of thin plate spline 
function and high interpolation accuracy by the incorporation of parametric linear sub-model, in 
addition to the independent spline variables; this software has been widely used in meteorological 
variables (Fick and Hijmans, 2017; Guo et al., 2020). 

Soil dataset, including soil types, soil organic matter content, soil sand and clay contents, and 
the percentage of soil organic carbon (SOC), was obtained from the China Soil Characteristics 
Dataset of the National Cryosphere Desert Data Center (http://www.ncdc.ac.cn) that was 
calculated from the 1:1,000,000 soil database of the second national soil census in China. The 
30-m digital elevation model (DEM) data were downloaded from the Geospatial Data Cloud 
(http://www.gscloud.cn). Soil dataset and DEM were used to calculate Zs. 

All above datasets were projected to the UTM-WGS84 coordinate system and resampled to 1 
km spatial resolution by the nearest neighbor resampling. Moreover, datasets for the whole year 
and in the growing season were calculated by the monthly average data of January through 
December and May through September, respectively; while for the non-growing season, datasets 
were calculated by the monthly average data of October through April of the next year. 

The 1-km gridded land use cover datasets in 2000 and 2020 were derived from the Resource 
and Environment Science and Data Center (http://www.resdc.cn). There were six categories of 
land cover types in this study: cropland, forestland, grassland, water body, bare land, and built-up 
land (John et al., 2009). 


2.3 Construction of modified remote sensing ecological index (MRSEI) 


Considering the regional characteristics of the dryland ecosystem in Otindag Sandy Land, we 
selected the vegetation status, water and heat conditions, and soil erosion to represent the local 
indicators of ecological quality in Otindag Sandy Land. The four indicators include Ze, Iw, In, and 
Is. Figure 2 shows the flowchart of spatiotemporal analysis of ecology quality in Otindag Sandy 
Land with MRSEI in this study. 

Vegetation cover is not only an important environmental element but also an indicator of the 
sensitive state of terrestrial ecosystems, and plays an important role in the analysis and evaluation 
of regional ecological and environmental factors (Tian and Min, 1998; Xu et al., 2019). J, reflects 
the growth status and distribution characteristics of green vegetation in the study area, and is 
expressed by Equation 1. 


I, = (NDVI-NDVI,9i1) f (NDVI jeg -NDVI i k (1) 


where NDVIsoii and NDVIveg are the NDVI values for the areas covered with bare soil and 
complete vegetation, respectively. 

Iw is expressed by the ratio of regional total precipitation to total potential evapotranspiration, 
which can effectively characterize the regional wet and dry conditions to determine the water 
requirements of ecosystems (Fernandes et al., 2016). It was calculated by the following equation: 


I, = PRE/ETo, (2) 
where, PRE is the total precipitation (mm); and ETo is the total potential evapotranspiration (mm), 


which was calculated based on the Penman-Monteith model using the interpolated climate 
datasets (Allen et al., 1998; Chen, 2021). 
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Fig. 2 Flowchart for constructing modified remote sensing ecological index (MRSEI) and spatiotemporal 
analysis of ecology quality with MRSEI in Otindag Sandy Land. MOD13A3, Monthly Normalized Difference 
Vegetation Index (NDVI) products from the Moderate-resolution Imaging Spectroradiometer (MODIS); DEM, 
digital elevation model; GEE, Google Earth Engine; ANUSPLIN, Australian National University Spline; NDVI, 
Normalized Difference Vegetation Index; ETo, total potential evapotranspiration; Jz, greenness index; Jw, wetness 
index; Zn, heatness index; Js, soil erosion index; PCA, principal component analysis; RSEI, remote sensing 
ecological index; ARSEI, arid remote sensing ecological index. 


In is closely related to the energy exchange processes within the ecosystem and reveals the 
potential amount of vegetation cover, which is expressed by the land surface temperature, an 
important indicator of the surface environment (Coutts et al., 2016). In this study, we replaced the 
land surface temperature in RSEI derived from remote sensing data (Xu et al., 2018) by air 
temperature from the National Meteorological Science Data Center (http://data.cma.cn). 

I; was calculated by the Revised Universal Soil Loss Equation (RUSLE) model, which 
synthetically considers the effects of precipitation and runoff on near-surface soil; this model has 
been widely used for the assessment of water and soil erosion (Toy and Osterkamp, 1995; 
Masroor et al., 2022). The indicator J; was calculated as follows: 


i= 4.4, (3) 
Ap = RXKXLXS, (4) 
A,=RxXKXLxSxCxP, (5) 
_wv4 n N 1.7265 
R =A Dex Pale )/n), (6) 
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where, Ap and A, are the potential and actual soil erosion (t/(km°-a)), respectively; R is the factor 
of precipitation erosion (MJ-mm/(km?-h-a)); K is the factor of soil erodibility (t-h/((MJ-mm)); L and 
S are the factors of the length and angle of the slope, respectively (dimensionless); C is the factor 
of vegetation cover and management; and P is the factor of soil and water conservation measures. 
Pa;; x represents the j” daily precipitation with erosiveness on the k half a month of the i” year, 
where k is the serial number of half a month in a year (A=1, 2, ..., 24) and i is the serial number of 
years (i=l, 2, ..., n). The j and N are the serial number and cumulative frequency of daily 
precipitation with erosiveness on the k half a month of the i” year, respectively. Kepic is the 
factor of soil erodibility in pre-revision. mc, msit, and ms represent the percentage contents of clay 
(<0.002 mm), silt (0.002—0.050 mm), sand (0.050—2.000 mm), respectively, with the unit of %. 
SOC is the percentage of soil organic carbon (%). 4 is the horizontal slope length (m), m is the 
slope length index, and @ is the slope degree (°). The elaborate calculation method of model 
factors refers to the study of Masroor et al. (2022). 

The above four indicators (le, Iw, In, and Zs) were projected to the UTM-WGS84 coordinate 
system at the spatial resolution of 1 km and temporal resolution of 1 month during 2001-2020, 
and then normalized to the dimensionless form in the range of 0-1. 

The weight for initial MRSEI (denoted as MRSEIo) based on the aforementioned four 
normalized indicators from 2001 to 2020 was calculated through PCA. This approach effectively 
avoids the bias of results caused by different artificial defined weights, making the results more 
objective and reliable (Hang et al., 2020). During 2001-2020, the contribution rates of the first 
component (PC1) were all greater than 73.00%, indicating that PC1 had integrated most of the 
characteristics of the four indicators (Fig. 3); therefore, PC1 was chosen to calculate the MRSEIo. 
The formula of MRSElb is expressed as follows: 


MRSEIy =1-| PCI, Ly ty4,) |, (12) 


g? w? 
where, MRSEIo is the initial value of MRSEI (from —1.000 to 1.000) at different time scales (i.e., 


the whole year, growing season, and non-growing season) for the given year. 
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Fig. 3 Contribution rates of the first principal component (PC1) of MRSEI indicators (greenness, wetness, 


heatness, and soil erosion) at the whole year, growing season, and non-growing season scales from 2001 through 
2020 


To facilitate the comparison among different time periods, we normalized the MRSEIo values 
to obtain the final MRSEI value based on Equation 13. 


MRSEI = (MRSEIy — MRSElmin )/(MRSElomax ~ MRSElpmin )> (13) 
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where MRSElomax and MRSElomin are the maximum and minimum values of MRSElbo, respectively. 
The higher MRSEI values indicate better ecological quality, while the lower MRSEI values 
represent poorer ecological quality. According to the Technical Criterion for Ecosystem Status 
Evaluation issued by the Ministry of Environmental Protection of China (Wang et al., 2021), we 
divided MRSEI (0.000—1.000) into five grades (worst: 0.000—0.200; poor: 0.200—0.400; general: 
0.400—0.600; good: 0.600—0.800; and best: 0.800—1.000) for further analysis. 


2.4 Statistical analysis 


2.4.1 Path analysis 


Path analysis is widely used to explore causal structure patterns between dependent and 
independent variables based on factor analysis and multiple linear regression (Fan et al., 2016; 
Shao et al., 2022). In this study, we regarded MRSEI as the dependent variable and Tg, Iw, In, and 
I; as the independent variables to clarify the correlations between them and initially verify the 
reliability and effectiveness of MRSEI. The larger the absolute value of path coefficient, the 
stronger the effect of the independent variable on MRSEI. The plus or minus path coefficient 
means its promotion or inhibition of MRSEI, respectively. 

2.4.2 Trend analysis 


Theil-Sen's slope method is a very useful nonparametric method for estimating linear trend (John 
et al., 2013; Xu et al., 2019). Therefore, it was used in this study to calculate the slope of each 
variable in each pixel and analyze the distribution pattern of the change trend (John et al., 2013). 
The slope (£) represents the change rate of variables. If the slope is negative (<0), it indicates 
that the variable is decreasing; otherwise, the variable is increasing ($>0). In addition, the 
non-parametric Mann-Kendall test (Niu et al., 2023) was used to evaluate the significance of the 
change trend of variables. The trend of MRSEI in Otindag Sandy Land was categorized into four 
types: significant increase (6>0 and P<0.05), increase (6>0 and P>0.05), significant decrease 
(B<0 and P<0.05), and decrease (6<0 and P>0.05). 

2.4.3 Barycenter model 


Because the dynamic variation of the barycenter can reflect the contrast and shift of regional 
element distribution (Huang et al., 2021), it was used to perform spatial analysis of barycenter 
changes in this study. The barycenter of MRSEI was calculated using the following equations: 


aaa x xCRR /(S) 5 R) (14) 


y=») mi 2 ERR, (X; _,CRR;] (15) 


where, X and Y are the coordinates of latitude and eek of b barycenter, respectively; x; and 
y; are the coordinates of latitude and longitude corresponding to the i and j" pixels, respectively; 
I and J are the number of pixels of latitude and longitude, respectively; CRR; and CRR; are the 
MRSEI values corresponding to the i and j® pixels, respectively. 

All data processing and analysis were executed in Python 3.7 (Anaconda Inc., Austin, USA), 
and the visualization implementation was conducted in OriginPro 2023 (OriginLab Inc., 
Northampton, USA) and ArcMap 10.2 (Esri Inc., Redlands, USA). 


2.5 Validation 


To further assess the feasibility of MRSEI, we compared MRSEI with the traditional RSEI (Xu, 
2013) and ARSEI for arid areas (Wang et al., 2020) based on 1-km MODIS datasets. RSEI was 
constructed with four indicators, including greenness, wetness, dryness, and heatness, whereas 
ARSEI took five indicators into account, comprising greenness, wetness, heatness, salinity, and 
land degradation. In this study, we calculated RSEI and ARSEI in the growing season of 2001, 
2005, 2010, 2015, and 2020 based on PCA via the Google Earth Engine. More information about 
the principles of RSEI and ARSEI can be found in Xu (2013) and Wang et al. (2020), 
respectively. 
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3 Results 


3.1 Feasibility analysis of MRSEI 


The four indicators (Ig, Iw, In, and Zs) of MRSEI were calculated by PCA on the basis of their 
contribution rates. Because the PC1 of the four indicators from 2001 to 2020 can integrate the 
most of the indicator characteristics, we took 2001 and 2020 as typical years with more detailed 
analysis. The subsequent path analysis further verified the significant correlation between the four 
indicators and MRSEI. 

As shown in Table 1, the PC1 loads of J, Iw, and J; were positive, and the average PC1 load of 
Iw (0.63) was larger than those of J, (0.06) and J, (0.35). This indicated that all these three 
indicators had positive effects on MRSEI, and the contribution rate of Iw was higher than those of 
I; and Iz, which was consistent with the results from path analysis. Furthermore, it was confirmed 
with the actual situation that precipitation is the key environmental factor in this area. Conversely, 
the PC1 loads of J, were negative, and the absolute value of its average load (0.68) was greater 
than those of Ie, Iw, and Zs. This indicated that J, had a negative effect on MRSEI, and its 
contribution to MRSEI was higher than Z, Iw, and Zs, also consisting with the results from path 
analysis. In addition, the absolute value of the PC1 load for J, was the highest in 2001, but the 
absolute value of the PC1 load for Zw was the highest in 2020, indicating that the main influencing 
indicator of MRSEI changed from /; in 2001 to Zw in 2020. The values of PC1 loads for the whole 
year were roughly close to those in the growing season, implying that ecological quality was 
affected more in the whole year than in the growing season. Therefore, it can be concluded that 
MRSEI based on PC1 is feasible and can be used to conduct the ecological quality assessment of 
Otindag Sandy Land. 

As shown in Figure 4 based on path analysis, Iw was the most important indicator affecting 
MRSEI with a highly significant correlation between them in all time periods from 2001 to 2020, 
followed by e, In, and J, in sequence. J, and MRSEI were negatively correlated, and the 
correlation was strong in the non-growing season. J, was positively and significantly correlated 
with MRSEI, with correlation coefficients of 0.80 in the growing season and 0.66 in the 
non-growing season, indicating that the influence of vegetation status (greenness) was weakened 
in the non-growing season. The results of path analysis also verified the feasibility of MRSEI 
from another perspective. 


Table 1 PC1 of MRSEI indicators (Zg, Iw, Jn, and Js) in 2001 and 2020 


Year Time scale os Eigenvalue Accumulated contribution 
L I, I, L rate (%) 

Whole year 0.32 0.54 -0.78 0.04 0.14 87.72 
2001 Growing season 0.40 0.57  —0.72 0.06 0.11 84.10 
Non-growing season 0.31 0.47 -0.82 0.06 0.06 74.46 
Whole year 0.37 0.73  —0.57 0.06 0.08 73.95 
2020 Growing season 0.44 0.73 -0.52 0.06 0.08 74.53 
Non-growing season 0.24 0.72 -0.65 0.06 0.06 73.96 


Note: PC1, first principal component of principal component analysis; MRSEI, modified remote sensing ecological index; J,, 
greenness index; /,,, wetness index; J», heatness index; /,, soil erosion index. 


3.2 Spatial changes of ecological quality 


The average MRSEI values in Otindag Sandy Land from 2001 to 2020 at different time scales are 
shown in Figures 5 and 6. The regions with the worst and poor grades of MRSEI at the whole 
year scale were scattered in the west of Otindag Sandy Land (Fig. 5al), accounting for 10.55% 
and 17.29% of the total area on average, respectively. As shown in Figure 5b1 and cl, the regions 
with these grades (worst and poor) shrunk in the growing season (8.92% and 15.96% of the total 
area, respectively) but expanded to the eastside of Otindag Sandy Land in the non-growing 
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Fig. 4 Path analysis of MRSEI and its four indicators (Zg, Iw, In, and Zs) in Otindag Sandy Land at the whole year 
(a), growing season (b), and non-growing season (c) scales from 2001 through 2020. Value around the double 
arrows indicates the correlation coefficient between independent variables, while value around the single arrow 
represents the path coefficient between independent and dependent variables. *, passing the 0.05 significant test; 
** passing the 0.01 significant test. 
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Fig. 5 Spatial distributions of multi-year average MRSEI and the corresponding MRSEI change trend at the 
whole year (al and a2), growing season (bl and b2), and non-growing season (cl and c2) scales from 2001 to 
2020 
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Fig. 6 Area proportions of different MRSEI grades in Otindag Sandy Land at the whole year (a), growing 
season (b), and non-growing season (c) scales from 2001 to 2020 


season (15.40% and 21.67% of the total area, respectively). In addition, the increased areas with the 
worst grade of MRSEI in the non-growing season expanded to most regions in Sonid Right Banner 
(Fig. 5c1). On the other hand, the regions with the good and best grades of MRSEI were mostly 
found in the southeastern part of Otindag Sandy Land, which accounted for 31.19% and 19.41% 
of the total area, respectively (Fig. 5b1), but the regions with the best grade of MRSEI were 
transformed to the good and general grades in the non-growing season (Fig. 5c1). This resulted in 
the regions with the best grade of MRSEI decreasing to 8.58% of the total area in the 
non-growing season. 

Correspondingly, the variations of MRSEI across Otindag Sandy Land varied considerably in 
space over the past 20 years. MRSEI at the whole year scale increased in 96.56% of the total area, 
with 75.62% of the total area showing a significant increasing trend (P<0.05; Fig. 5a2). The 
increase of MRSEI in the western and southern parts of Otindag Sandy Land mainly existed in the 
growing season (Fig. 5b2), while the increase of MRSEI in the northern part and eastside of 
Otindag Sandy Land mostly existed in the non-growing season (Fig. 5c2). In contrast, MRSEI at 
the whole year scale decreased in the remaining regions (3.44% of the total area), mostly 
distributed in Taibus Banner of the southern part of Otindag Sandy Land, mainly due to its 
ecological degradation in the non-growing season. In addition, the regions with the decreased 
MRSEI reached 18.64% of the total area in the non-growing season, which were scattered in the 
southern, western, and northern parts of Otindag Sandy Land; however, only 0.27% of the total 
area showed significant decrease in MRSEI (P<0.05). The comparison of MRSEI in 2001, 2005, 
2010, 2015, and 2020 further indicated that ecological quality in the growing season was stably 
improved in the whole Otindag Sandy Land, while MRSEI in the non-growing season showed a 
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stronger inter-annual fluctuation (Fig. S1). Severe ecological quality degradation was found in the 
non-growing season of 2010, which was mainly occurred in the southern and northern parts of 
Otindag Sandy Land (Fig. S1). Overall, it appeared that ecological quality over most regions of 
Otindag Sandy Land had improved from 2001 to 2020, but there were still some small patches 
showed degradation in ecological quality. 

Table 2 exhibits the area transfer matrix of MRSEI grades at the whole year scale from 2001 to 
2020. Relative to 2001, the area with the best, good, and general grades of MRSEI increased 
while the area with the worst and poor grades of MRSEI decreased in 2020. The area with the 
best and good grades of MRSEI increased by 20,493.89 km? (from 30,214.83 km? in 2001 to 
47,105.73 km? in 2020 for the best grade and from 42,136.76 km? in 2001 to 45,739.75 km? in 
2020 for the good grade, accounting for 9.77% and 2.09% of the total area, respectively), with 
most converted from the area with the general grade of MRSEI (19,351.89 km?), and 1142.00 km? 
from the area with the worse and poor grades of MRSEI. Meanwhile, the area degrading from the 
good grade of MRSEI to the general grade was 2254.99 km’. The area with the general grade of 
MRSEI in 2020 increased by 10,668.94 km?, accounting for 6.17% of the total area. The area with 
the worst grade of MRSEI was significantly decreased by 27,357.84 km? (from 35,832.79 km? in 
2001 to 8474.95 km? in 2020), indicating that only 4.90% of Otindag Sandy Land had the worst 
grade of MRSEI in 2020. 


Table 2 Area transfer matrix of MRSEI with different grades in Otindag Sandy Land from 2001 to 2020 (unit: 


km?) 
Worst Poor General Good Best Total (2001) 
Worst 8462.95 18,262.89 8648.95 458.00 0.00 35,832.79 
Poor 12.00 3571.98 21,380.88 684.00 0.00 25,648.86 
General 0.00 9.00 17,351.23 21,541.88 65.00 38,967.11 
Good 0.00 0.00 2254.99 22,205.87 17,675.90 42,136.76 
Best 0.00 0.00 0.00 850.00 29,364.83 30,214.83 
Total (2020) 8474.95 21,843.87 49,636.05 45,739.75 47,105.73 172,800.35 


3.3 Temporal evolution of ecological quality 


Figure 7a shows the temporal variations in the ecological quality of Otindag Sandy Land from 
2001 to 2020. MRSEI in the growing season exhibited a significant increasing trend with the 
change rate of 0.004/a (P<0.05) from 2001 to 2020. The corresponding trend in the non-growing 
season slowed down after 2010 with insignificant change rate (P>0.05). In addition, the 
inter-annual fluctuation was high in the non-growing season, such as the dramatic decline from 
2015 to 2018. Due to the dominance of MRSEI in the growing season, the change of MRSEI at 
the whole year scale followed the change trend of MRSEI in the growing season; specifically, it 
increased significantly with the change rate of 0.006/a (P<0.05) from 2001 to 2020. 

Figure 7b—d depicts the barycenter movements of MRSEI from 2001 to 2020. The longitudinal 
barycenter of MRSEI at the whole year scale declined (P<0.05) and tended to move by 0.2° 
toward the west. However, the latitudinal barycenter was barely unchanged and did not pass the 
significant test (P>0.05). It can be seen that the longitudinal and latitudinal barycenter changes of 
MRSEI at the whole year scale were basically determined by the shifts in the growing season. The 
longitudinal and latitudinal barycenter of MRSEI in the non-growing season showed an 
insignificant increasing trend (P>0.05), with the longitudinal barycenter inclining back to the east 
and the latitudinal barycenter at a higher latitude. The longitudinal movement direction of the 
barycenter in the non-growing season was opposite to that in the growing season in some years 
(e.g., 2005 and 2015). On the whole, the barycenter of MRSEI in Otindag Sandy Land tended to 
move westward, which indicated that in the past 20 years, ecological quality in the western part of 
Otindag Sandy Land has improved significantly. 
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Fig. 7 Temporal change of MRSEI in Otindag Sandy Land (a) and its barycenter movements (b—d) at the whole 
year, growing season, and non-growing season scales from 2001 to 2020 


3.4 Impact of climate and land cover on ecological quality 


Figure 8al and a2 depicts the relationships of MRSEI changes with the variations in precipitation 
and air temperature. Precipitation was positively related to MRSEI; specifically, precipitation 
showed an increment of 1.5 mm/a, corresponding to an increase of MRSEI by 0.002 (+0.002) per 
year (Fig. 8al). In contrast, air temperature was negatively related to MRSEI (Fig. 8b1); 
specifically, air temperature decreased by 0.02°C/a, corresponding to an increase of MRSEI by 
0.002 (40.003) per year. However, precipitation and air temperature were not linearly related to the 
changes in MRSEI. The degree of MRSEI changes became more apparent when the increase of 
precipitation was higher than 4.5 mm/a and the decrease of air temperature was more than 
0.02°C/a. As seen from Figure 8a2, precipitation showed increasing trends over the whole Otindag 
Sandy Land, with the change rates of 1.5-4.5 mm/a in more than 76.49% of the total area, whereas 
the increment of precipitation of more than 4.5 mm/a was only found in the middle of Sonid Right 
Banner and Sonid Left Banner in the western part of Otindag Sandy Land. Correspondingly, air 
temperature decreased in the same regions where the increment of precipitation was more than 4.5 
mm/a (Fig. 8b2). This synergistic effect directly increased MRSEI, which significantly improved 
ecological quality from the worst grade to poor or general grade (Figs. 5 and S1). 

Land cover types had different effects on MRSEI. For the four major land cover types in 
Otindag Sandy Land (forestland, cropland, grassland, and bare land), the multi-year average 
MRSEI was 0.805, 0.779, 0.544, and 0.433, respectively (Fig. 8c1). For land cover changes (Fig. 
8c2), MRSEI increased by 0.004 (40.004) per year for both the conversion of cropland to 
grassland and the inversion of grassland to cropland. The transformation of bare land to grassland 
increased MRSEI by 0.006 (+0.004) per year. Afforestation was the most effective forcing to 
improve MRSEI in the study area, with MRSEI increased by 0.009 (+0.004) per year. 
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Fig. 8 Relationships of the slope of MRSEI with the slopes of precipitation (al) and air temperature (b1), and 
spatial distributions of changes in precipitation (a2) and air temperature (b2), as well as multi-year average 
MRSEI over different land cover types (cl) and MRSEI changes with transformation of land cover types (c2). 
Bars mean standard errors; dot means the average of the slope of MRSEI. Colors in Figure 8al and b1 correspond 
to color classes in Figure 8a2 and b2, respectively. 


4 Discussion 


4.1 Validation of the effectiveness of MRSEI 


We examined the ability of MRSEI to capture the local texture information of water bodies and 
topography (Fig. 9). Figure 9a showed that MRSEI could accurately capture the shape of Dalai 
Nur Lake and the presence of some small lakes, as well as the higher vegetation coverage in the 
east of Otindag Sandy Land, where MRSEI values were also higher. Meanwhile, MRSEI 
accurately captured the annular vegetation coverage zone with the change of topography, as well 
as the buildings in the southeastern part of annular zone that corresponded to lower MRSEI values 
(Fig. 9b). Clearly, MRSEI could well reflect the actual scene of land surface in Otindag Sandy 
Land. 

Comparisons among the contribution rates of PC1 of RSEI, MRSEI, and ARSEI indicators in 
the growing seasons of 2001, 2005, 2010, 2015, and 2020 are shown in Figure 10a. The 
contribution rates of PC1 of MRSEI indicators were the highest in all years, followed by RSEI 
and ARSEI indicators. Furthermore, the spatiotemporal variations were compared among these 
three indices (Fig. 10b). RSEI, ARSEI, and MRSEI depicted similar low-high zonal changes, with 
values increasing from the west to the east and culminating at nearly 117°-118°E, but they 
showed different spatial heterogeneity, with MRSEI being the highest and RSEI being the lowest. 
Interestingly, all these three indices exhibited increasing trends at the temporal scales, but their 
change rates were slightly different. Overall, all these three indices reflected that ecological 
quality was improved in Otindag Sandy Land over the past 20 years. However, RSEI did not 
consider the local characteristics of Otindag Sandy Land (e.g., soil erosion) in drylands even 
though the contribution rates of PCl of RSEI indicators were higher than 70.00%. The 
contribution rates of PC1 of ARSEI indicators were between 56.50%—68.38%. In this case, if only 
the PC1 was considered, much information may be omitted (Jiang et al., 2021). Therefore, 
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MRSEI can be considered as an effective alternative index for ecological quality assessment in 
Otindag Sandy Land. Our results also showed that the selection of indicators should ideally be 
adapted to the local contradictions so that it can accurately characterize the changes in regional 
ecological quality. 


(a) Study area N (b) TCC and MRSEI in 2020 
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— : = MRSEI 
c 1.000 
0 100km 
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Fig. 9 Local validation of MRSEI in capturing the texture information of water bodies (a and b) and topography 
(a and c) between MRSEI and real scenes of surface in 2020. TCC is the true color composite image (RGB 432) 
of Landsat 8 filtered; it was calculated by the median composite method of the images with the cloud cover of less 
than 5% in the growing season (from May to September) in 2020 by the Google Earth Engine. 
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Fig. 10 Comparison between the contribution rates of PC1 of RSEI, MRSEI, and ARSEI indicators in the 
growing seasons of 2001, 2005, 2010, 2015, and 2020 (a) and zonal distribution and temporal changes (the inset) 
of RSEI, MRSEI, and ARSEI (b). *, passing the 0.05 significant test. Bars mean standard errors and dot means the 
average of slope. 


4.2 Impact of climate and land cover on MRSEI 


Based on the results of PCA and path analysis, the wetness indicator played a major role in 
improving ecological quality, which confirmed the pivotal status of water resources in arid and 
semi-arid areas (Wu et al., 2022). In Otindag Sandy Land, precipitation in the eastern and 
northeastern parts is higher (Ma et al., 2017), which is consistent with the spatial distribution of 
MRSEI. The western part of Otindag Sandy Land with less precipitation and sparser vegetation 
(Shi et al., 2004) was rated as the region with the worst grade of MRSEI. The fragile ecosystems 
in Otindag Sandy Land are sensitive to water resources, and spatial differences in water resources 
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can lead to the heterogeneity of ecological quality in Otindag Sandy Land. 

The overall ecological quality in Otindag Sandy Land was improved from 2001 to 2020, but it 
showed large inter-annual variations. The improved regions accounted for approximately 96.56% 
of the total area across Otindag Sandy Land. The main reason is that climate in Otindag Sandy 
Land becomes more humid (Xu et al., 2020; Zhao et al., 2021), as shown in Figure 8, which is 
conducive to relieving local water resources shortage. Increased precipitation tends to improve 
vegetation coverage, especially in the western part of Otindag Sandy Land with grasslands and 
bare lands (Cong et al., 2022), where MRSEI increased significantly and ecological quality was 
improved (Fig. 5). However, the reported warming in the non-growing season might cause the 
evapotranspiration increased and soil dried out (Zhang et al., 2019; Xu et al., 2020). The 
improvement process of vegetation and soil conditions in the growing season were interrupted or 
deteriorated in the following non-growing season (Fig. 5). These results indicated that ecological 
quality has seasonal changes, and attention should be paid to its seasonal degradation. Given that 
future climate will get more humid in this region, there is still room for sustainable ecological 
improvement in Otindag Sandy Land from the water supply perspective. 

In addition to the influence of climate change, the implementation of ecological restoration 
measures in Otindag Sandy Land, such as maintaining grass and livestock balance, returning 
farmland to forests, and implementing enclosures and grazing prohibitions (Yuan et al., 2016; Ma 
et al., 2017), can be conducive to improving the ecological quality of Otindag Sandy Land. Due to 
the effectiveness of these measures, regional vegetation has been effectively restored, which 
alleviated ecological deterioration to some extent, especially for the grasslands and bare lands. 
Combined with the increase in precipitation, many regions with poorer ecological quality (worst 
and poor grades) were transformed into regions with better ecological quality (general, good, and 
best grades). Afforestation in the forest regions of the eastern part of Otindag Sandy Land with 
richer precipitation directly increased NDVI, in turn, enhanced ecological quality therein. 
However, afforestation may lower the groundwater level in arid areas (Chen et al., 2018; Lu et al., 
2018). We therefore recommend that local climate change and water resources should be taken 
fully into account when implementing ecological restoration measures. 

In terms of land cover and land use change, grassland is the dominant land cover type in 
Otindag Sandy Land; the quality of grasslands has been improved over the last 20 years due to 
increased precipitation and ecological restoration measure implementations. The land cover type 
with denser vegetation had higher MRSEI. However, this does not mean that the changes of 
higher MRSEI in the regions with denser vegetation are positive. For example, higher MRSEI in 
the southern part of Otindag Sandy Land decreased, especially in the non-growing season. As 
shown in Figure 1, the transformation of grasslands to croplands as well as the expansion of 
croplands that occurred in the southern part of Otindag Sandy Land may lead to soil salinization 
and degradation of natural habitats (Tang et al., 2021; Zhuang et al., 2021), which can explain the 
degradation in the southern part of Otindag Sandy Land with good ecological quality, especially 
in Taibus Banner. In sum, ecological restoration is important for sustainable development, but 
local economic benefits from agriculture and husbandry expansion must not be ignored. Thus, 
how to balance local economic benefits and ecological revitalization is the most crucial factor for 
regional ecological sustainability in the future. 

Although the MRSEI changes in Otindag Sandy Land showed that ecological quality is 
gradually improving, there are still uncertainties in restoring ecological quality in the future. 
According to previous studies (e.g., Chen et al., 2020), these highly precipitation-sensitive 
regions are ecologically fragile and easily suffer from droughts because of rain-fed vegetation 
properties, so both the significant inter-annual and seasonal fluctuations of precipitation and the 
increased water consumption due to rapid warming can increase the threat of short-term drought 
to its local ecosystems (John et al., 2013). Thus, the seasonal evaluation of ecological conditions 
as well as corresponding conservations are very important. How to reasonably and scientifically 
formulate measures to improve the ecological quality of Otindag according to local characteristics 
and environmental conditions needs further evaluations. From this study, the priority is to 
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maintain the natural landscapes and avoid overgrazing in the western part of Otindag Sandy Land 
with scarce water resources. In regions with severe ecological degradation, it is essential to 
implement graded governance and encourage local tourism to satisfy the harmonious 
development of economy and nature. On the other hand, since the development of MRSEI aims to 
explore the ecological quality of Otindag Sandy Land with typical dryland ecosystem, its 
practicality is temporarily limited to dryland ecosystems, and whether it is appropriate to extend 
to a broader area needs further research. 


5 Conclusions 


In this study, we constructed the MRSEI integrating greenness, wetness, heatness, and soil 
erosion indicators coupled with the local characteristics and environmental conditions based on 
PCA and path analysis to evaluate the ecological quality of Otindag Sandy Land from 2001 to 
2020. All contribution rates of PC1 were above 73.00% during 2001-2020. MRSEI based on PC1 
appeared promising in Otindag Sandy Land. Results of PCA and path analysis showed that 
greenness, wetness, and soil erosion indicators had positive effects on MRSEI, whereas heatness 
had a negative effect. The contribution rate of PC1 of wetness was the highest, followed by 
heatness, greenness, and soil erosion. 

In the past 20 years (2001-2020), the changes of ecological quality in Otindag Sandy Land 
showed seasonal and regional differences. Ecological quality in the growing season was improved 
significantly, but in the non-growing season it fluctuated with an insignificant increasing trend. 
Spatially, the regions with the good and best grades of MRSEI were mostly distributed in the 
eastern part of Otindag Sandy Land, while the regions with the worst and poor grades of MRSEI 
were concentrated in the western part of Otindag Sandy Land. However, the regions with the 
worst and poor grades of MRSEI respectively accounted for 8.92% and 15.96% of the total area 
in the growing season but 15.40% and 21.67% of the total area in the non-growing season. 
Ecological equality was improved in most regions of Otindag Sandy Land in the growing season 
but degraded in the western, northern and southern parts of Otindag Sandy Land in the 
non-growing season (18.64% of the total area). Altogether, the area with the worst grade of 
MRSEI shrunk by 15.83% of the total area during the study period, from 35,832.79 km? in 2001 
to 8474.95 km? in 2020, while the area with the best grade of MRSEI in the east increased by 
9.77% of the total area from 30,214.83 km? in 2001 to 47,105.73 km? in 2020. Over time, 
spatiotemporal variations of ecological equality varied with seasons. Such seasonal differences 
should be taken into account for better reflecting ecological quality assessment results. For 
multi-year analysis, we suggest to choose indicators to represent seasonal ecological conditions 
according to the local characteristics of the study area. 
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Fig. S1 Spatial distributions of MRSEI at the whole year (al—a5), growing season (bl-b5), and non-growing 
season (cl—c5) scales in 2001, 2005, 2010, 2015, and 2020 in Otindag Sandy Land 


